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The generation, propagation and radiation of sound from a perfectly expanded Mach 2.5 
cold supersonic jet flowing through an enclosed rigid-walled duct with an upstream J- 
deflector have been numerically simulated with the aid of OVERFLOW Navier-Stokes CFD 
code. A one-equation turbulence model is considered. While the near-field sound sources are 
computed by the CFD code, the far-fleld sound is evaluated by Kirchboff surface integral 
formulation. Predictions of the farfield directivity of the OASPL (Overall Sound Pressure 
Level) agree satisfactorily with the experimental data previously reported by the author. 
Calculations also suggest that there is significant entrainment of air into the duct, with the 
mass flow rate of entrained air being about three times the jet exit mass flow rate. 


Nomenclature 

c - sound velocity 

dj = jet exit diameter 

/ = frequency 

m = mass flow rate 

M = Mach number, ulc 

p = pressure 

p' = acoustic pressure disturbance 

r = distance from the sound source 

Re - Reynolds number, PjUjdj / pj 

St ~ pressure 

T = temperature 

u - velocity 

x ,y,z = Cartesian coordinates 

GREEK SYMBOLS 

p - dynamic viscosity 

p = density 

SUBSCRIPTS 

j = jet 

oo = ambient fluid 


I. Introduction 

1' he overall process of sound emission from partially ducted rocket exhausts is of great practical interest in the 
understanding of sound suppression systems for launch vehicles. Clean (dry) launch pads with ducted exhausts 
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are preferable to those fitted with water deluge systems for sound suppression with regard to maintenance costs and 
the frequency of launches. A detailed knowledge of the behavior of the sound radiation from these ducted exhausts 
is useful in the design and optimization of the sound suppression systems. 

The totality of the far field sound from ducted exhausts is composed of the transmission of sound from the free 
supersonic jet (upstream of the duct inlet) to the far field, transmission of sound within the duct, and radiation of 
sound from the duct exit and the subsonic jet (exiting the duct) to the far field. Recent experiments on enclosed duct 
by Kandula et al. 1 ’ 2 3 have indicated that the jet confinement has a significant effect on sound emission and directivity, 
suggesting that the duct is modifying the sound generation and propagation. According to Lighthill 4 ' 5 , turbulent 
mixin g is responsible for sound generation in jets, with Mach wave radiation (Fflowcs Williams; Tam) 
manifesting itself in supersonic turbulent jets due to supersonic convection of large scale eddies relative to the 
ambient. 

A number of studies were reported on enclosed subsonic flows typical of turbofan engine inlet/exhaust ducts 
(for example, Eversman & Roy 8 ; Ozyoruk & Long 9 ). However, many of the theoretical studies on sound emission 
from turbulent supersonic jets are focused on unconfined (free) jets (Freund et al. 10 ; Garnet & Estivalezes 11 ; Kandula 
& Caimi 12 ). The present study is concerned with a theoretical simulation of the effect of the enclosing duct on the 
radiated noise from cold supersonic jets. Predictions of sound power levels will be compared to the experimental 
data 1 ' 3 . 


I. Computational Simulation 


A. Geometry and Experimental Data 

Fig. 1 shows the geometry of the jet-duct configuration considered in the experimental data of Refs. 1-3. The 
origin of the coordinate system lies in the ground plane, with the y-axis coinciding with the jet axis. The jet 

impinges vertically on an upstream J-deflector (30 deg. to the vertical) of the duct before exiting horizontally 
through the duct. The nozzle exit diameter is 1 in. The nozzle exit is 73 in. above the ground level, and 10 in. above 
the duct inlet plane. The duct exit cross section is 6 in. x 12 in. Photographic views of the free jet and the ducted jet 
testing with Mach 2.5 cold nitrogen jet are displayed in Figs. 2a and Fig. 2b respectively. Fig. 2c shows the sound 
directivity in the ducted jet system 1 ’ 3 , with data taken along an 80 in. arc radius about the jet axis, and in a horizontal 
plane 54 in. above the ground level. In Fig. 2c, the quantity h refers to the height of the nozzle exit plane above the 
duct inlet plane. 

B. Grid System 

A front view of the grid system is presented in Fig. 3a. Fig. 3b shows a representative view of the lateral extent 
of the grid system. The overall grid system consists of 8 blocks with a total of about 4x1 0 6 grid points. Individual 
grids for the jet, duct, etc., are generated using GRIDGEN grid generation program. The inter-grid communication is 
established by the Pegasus Code (Benek et al. 13 ). 

C. Flow Parameters 

The Mach number of the exhaust jet is 2.5, and the jet Reynolds number is 58,300. The jet is perfectly expanded 
(nozzle exit pressure equals the ambient pressure of 14.7 psia). The exit static temperature of the jet is 222 R, 
indicating that it is a cold jet. The ambient temperature is 530 R. 

D. Solution Procedure 

The nearfield solution for the mean flowfield and the acoustic field is obtained by the OVERFLOW CFD 
Navier-Stokes code (Kandula & Buning 14 ). A one-equation turbulence model of Spalart-Allmaras 15 is employed, 
which solves for the turbulent kinetic energy. Multi-gridding scheme is chosen for accelerated convergence. 
Initially the steady hydrodynamic flowfield is obtained by running the code with time-stepping. Steady state is 
achieved in about 2000 iterations. Starting from the steady state solution, a three-dimensional periodic disturbance is 
applied at the nozzle exit. The disturbance frequency is based on a Strouhal number St =0.2. Thomson’s outflow 
acoustic boundary condition (Thompson 16 ) and radiation boundary condition of Tam and Webb 17 are considered. 
The code is run for about 22000 iterations to ensure that a periodic state is approached. 

In the far field, the sound pressure levels are determined by Kirchhoff formulation (Kirchhoff 18 with data on the 
Kirchhoff surface obtained by the CFD code. The method is presented in Lyrintzis & Mankbadi 19 , and has been 
recently considered by the author for an axisymmetric supersonic jet simulation (Kandula and Caimi 12 ). The data on 
the Kirchhoff surface includes time history of pressure, normal pressure derivative, and pressure-time derivative. 
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n. Results 


A. CFD/Kirchhoff Solution 

1. Air Entrainment 

Calculations show that there is significant entrainment of air into the duct. Fig. 4 displays the convergence 
history of the duct exit mass flow rate until steady state is established. It is seen that the entrained air flow is about 
3.2 times the mass flow rate of jet exit. This result is representative of entrainment in practical exhaust duct 
configurations. 

2. Instananeous Flow field 

Fig. 5 displays the instantaneous mach number contours in the duct and in the region downstream of the duct exit 
at a lateral location near the axis of the duct. While the jet exit Mach number is 2.5, the Mach number of flow at the 
duct exit is subsonic around Mach 0.45. 

3. Acoustic Presssure 

The acoustic pressure (obtained by subtracting the mean pressure from the instantaneous pressure) for the jet- 
duct system is illustrated in Fig. 6a. The directivity of Mach wave radiation from the initially supersonic jet 
(unenclosed portion) is observed. The plot also shows the diffraction pattern of sound around the duct. Fig. 6b shows 
an expanded view of the acoustic pressures within the duct, showing the character of sound reflections and 
propagation within the duct. It also reveals the acoustic pressure contours in the jet, showing the three-dimensional 
character of the disturbance in the jet. Calculations show that the OASPL at the location of the jet impingement 
within the duct is of the order of about 180 dB. Fig. 6c exhibits the acoustic pressure in a near lateral view, showing 
the near-spherical spreading of the subsonic exhaust from the duct. 

B. Comparison with Experimental Data 

The acoustic far field computations of sound levels are based on Kirchhoff surface integral. Fig. 7 presents a 
comparison of the predictions of the directivity of the OASPL with the test data 1 " 3 for a closed duct. The OASPL 
data are presented as a function of the azimuthal angle measured from the axis normal to the duct exit plane. While 
the free jet data for the OASPL shows symmetry, the OAPSL for the closed duct shows considerable directivity, 
with the OASPL decreasing with the increasing angle. Satisfactory agreement is noted between the theory and the 
data. The theory predicts the OASPL with an accuracy of about 2 dB. The results and comparisons will be useful in 
the development of scaling laws for ducted jet noise, such as those considered for free jets (Kandula and Vu 20 ). 

HI. Conclusion 

Numerical simulation of sound from a supersonic jet passing through a rigid duct with a J-deflector has been 
carried out. The results provide details of the sound radiation from the free jet and the acoustic field within the duct. 
The presence of the duct significantly modifies the character of the far-field directivity. Significant level of 
entrainment of air into the duct has been noted. 
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Fig. 1 Geometry of the ducted exhaust jet configuration. 



Fig. 2a. Photograph of the nitrogen supersonic free jet. 
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Fig. 2b. Photograph of the nitrogen supersonic jet flowing 
through an enclosed duct. 
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Mach 2.5 Jet 

Nozzle exit Diameter: 1 in. 

Nozzle Exit Plane: 73 in. above GL 
Acoustic Sensors: 80 in. arc radius; 54 in above GI 


-Free Jet (Run# 2) 

Run# 1 (Closed Duct, h- 10 in.) 
Run# 5 (Closed Duct, h« 5 in.) 
Run# 9 (Closed Duct, h*-l in.) 


Fig. 2c Effect of nozzle height on the directivity of sound 
radiation for a closed duct. 
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Fig. 3a Grid system for the exhaust (side view). 



Fig. 3b Grid system for the exhaust (lateral view). 
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Fig. 4 Convergence history of the duct exit mass flow rate. 




Fig. 5 Instantaneous Mach number contours in the exhaust jet. 
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Fig. 6a Acoustic pressure contours in the jet-duct system (excluding the duct). 



Fig. 6b Acoustic pressure contours in the duct region. 
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Fig. 6c Acoustic pressure in the lateral view (including the duct). 
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Fig. 7 Comparison of the directivity of Overall Sound Pressure Level. 
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